Hybrid model calculations of direct photons in high-energy nuclear collisions 
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Direct photon emission in heavy-ion colhsions is calculated within a relativistic micro+macro 
hybrid model and compared to the microscopic transport model UrQMD. In the hybrid approach, the 
high-density part of the evolution is replaced by an ideal 3-dimensional hydrodynamic calculation. 
This allows to examine the effects of viscosity and full local thermalization, in comparison of the 
transport model to the ideal fluid-dynamics. We study the origin of high-px photons as well as the 
impact of elementary high-y^ collisions. We further explore the contribution of different production 
channels and non-thermal radiation to the spectrum of direct photons. Detailed comparison to the 
measurements by the WA98-collaboration are also undertaken. 



o 

(N 



Oh' 

I ■ 

Oh- 

<D : 



m ■ 
> ■ 

\o ■ 

vn ■ 
o ; 
o\ . 
o . 



X 



I. INTRODUCTION 

Creating and studying high-density and -temperature 
nuclear matter is the major goal of heavy-ion experi- 
ments. A state of quasi-free partonic degrees of freedom, 
the Quark-Gluon-Plasma (QGP) [l|,0 may be formed, if 
the energy density in the reaction is high enough. Strong 
jet quenching, high elliptic flow and other observations 
made at the Relativistic Heavy Ion Collider (BNL-RHIC) 
suggest the successful creation of a strongly coupled QGP 
(sQGP) at these energies 0-0] • Possible evidence for the 
creation of this new state of matter has also been put for- 
ward by collaborations at the Super Proton Synchrotron 
(CERN-SPS), as for instance the step in the mean trans- 
verse mass excitation function of protons, kaons and pi- 
ous and the enhanced /v +/7r+-ratio Q. 

Out of the many possible observables, electromagnetic 
probes have the advantage of leaving the hot and dense 
region undisturbed: once they are created, they escape 
freely from the reaction zone, due to their negligible 
rescattcring cross-sections. Besides dileptons, direct pho- 
ton emission is therefore of greatest interest to gain in- 
sight into the early, hot and therefore possibly partonic 
stages of the reaction. Direct photons arc distinguished 
from the bulk of photons as those coming from collisions 
and not decays. 

Unfortunately, the overwhelming amount of photons 
in heavy-ion collisions comes from hadronic decays in the 
late stages, mostly tt^ — >■ 77. These decay-photons im- 
pose a serious challenge for the experimental extraction of 
direct photon data. Up to now, several experiments have 
gone through the challenge to obtain the spectra of direct 
photons: Helios, WA80 and CERES (all at CERN-SPS) 
could publish upper limits, while WA98 (CERN-SPS) Q 
and PHENIX (BNL-RHIC) have published explicit 

data points for direct photons. 

On the theoretical side, calculations for the elementary 
photon production processes are known since long, see 
e.g. Kapusta et al. 11 [ and Xiong et al. [l^ . The major 
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problem here is the difficulty to describe the time evolu- 
tion of the produced matter, which is, up to now, not 
possible from first principle Quantum Chromodynam- 
ics (QCD). One has to rely on well-developed dynamical 
models to describe the space-time evolution of the nuclear 
interactions in the hot and dense stage of the reaction. 
A well-established approach to explore the dynamics of 
heavy-ion reactions is relativistic transport theory [isl - 
[2lj | . In this kind of microscopic description, the hadronic 
and/or partonic stage of the collision is described under 
certain approximations. Most transport models, for in- 
stance, cannot describe collisions with more than two in- 
coming particles, which restricts the applicability to low 
particle densities, where multi-particle interactions are 
less important. Some attempts to include multi-particle 
interactions do exist [H, [22| - [26j , but this field of study is 
still rather new. The coupling of a partonic phase with 
a hadronic phase poses another challenge on transport 
models, because the microscopic details of that transi- 
tion are not well known. Another complication in the 
transport approach is that all microscopic scatterings are 
explicitly treated in the model and therefore the cross- 
sections for all processes must be known or extrapolated. 
However, for many processes high quality experimental 
data are not available, and therefore a large fraction of 
the cross-sections have to be calculated or parametrized 
by additional models. 

Relativistic, (non-)viscous fluid- or hydrodynamics is a 
different approach to explore the space-time evolution of 
a heavy- ion collision ^ ^] . It constitutes a macroscopic 
description of the matter that is created, assuming that 
at every time and in every place the matter is in local 
thermal equilibrium. This assumption can only be true 
if the matter is sufficiently dense, so in the late stages of 
a heavy-ion collision fluid-dynamics looses applicability. 
In addition, the requirement of local thermal equilibrium 
restricts the starting time of the hydrodynamic model. 
An advantage is that in the dense stages, hydrodynamics 
can propagate any kind of matter, and also allows for 
transitions between two types of matter, e.g. QGP and 
hadron gas, if an appropriate Equation of State (EoS) is 
provided. Fluid-dynamics can therefore be used to study 
hadronic and partonic matter in one common framework. 
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The restrictions of these kinds of models can be loos- 
ened as well. By introducing viscosity and heat conduc- 
tivity, perfect thermal equilibrium does not have to be 
present at any point. However, even with second order 
corrections the matter has to be close to equilibrium [43l - 

Input to solve the hydrodynamic differential equations 
are the boundaries, i.e. the initial state (the distributions 
of all relevant densities and currents at the time the evo- 
lution starts), the Equation of State providing the pres- 
sure as function of the energy and baryon densities, that 
describes the behaviour of the matter that is considered, 
and the freeze-out hypersurface. 

Finally, approaches to the theoretical description of 
direct photon spectra may include calculations in per- 
turbative Quantum-Chromo-Dynamics (pQCD). Calcu- 
lations based on pQCD describe the high-p^ photon data 
in proton-proton collisions very well and, if scaled by the 
number of binary nuclcon-nucleon collisions, also those in 
heavy-ion reactions. The range of applicability of these 
calculations is, however, limited to high transverse mo- 
menta ^ 1 GcV. 
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FIG. 1: (Color Online.) UrQMD-cross-sections for tt+tt"- 
collisions as function of center-of-mass energy. We show the 
resonant hadronic cross-section (red dashed line), the cross- 
section for the formation of strings (blue dotted line) and for 
hard scatterings via PYTHIA (green dash-dotted line). The 
peak at the p-meson pole mass has been cut out for better 
visibility. 



II. THE HYBRID MODEL 



A. Transport model 



The different approaches have different applicability 
limits. E.g., thermal rates can only be applied if the as- 
sumption of local thermal equilibrium is fulfilled. Photon 
emission rates can then be calculated by folding the par- 
ticle distribution functions of the participating particle 
species with the respective cross-sections. This frame- 
work can be applied to either static models, simplified 
hydrodynamics-inspired models such as the blast wave 
model and to full fluid-dynamic calculations. The space- 
time evolution of a reaction as predicted by microscopic 
theories can be averaged over in order to app ly thermal 
rates to the coarse-grained distributions |39|. The ap- 
plication of microscopic cross-sections can only be un- 
dertaken in a model where all microscopic collisions are 
known. That limits the field of use to transport models. 
For previous calculations of photon spectra from trans- 
port models see e.g. jiol - lisj . 

In this paper, we investigate the spectra of direct 
photons coming from microscopic hadronic scatterings, 
thermal hadronic and partonic emission and hard ini- 
tial pQCD scatterings. We compare results from a 
purely microscopic model to those from an integrated 
micro-f-macro hydrodynamic approach that embeds a hy- 
drodynamic phase into the UrQMD approach. The pa- 
per is organized as follows. First, we explain the hybrid 
model (Section |ll|. In the following Section HTll we elab- 
orate on the photon sources considered in our model. In 
Sections |lVl |V] and |Vll we apply our model to compare 
thermal rates from microscopic theory to those from the 
literature, compare different physics assumptions with 
experimental data from the WA98 collaboration Q and 
analyse the sources of photon emission. 



UrQMD v2.3 (Ultra-relativistic Quantum Molecular 
Dynamics) is a microscopic transport model [l3, E, . 
It includes all hadrons and resonances up to masses 
TO 2.2 GeV and at high energies can excite and frag- 
ment strings. The cross-sections are either parametrized, 
calculated via detailed balance or taken from the additive 
quark model ( AQM) , if no experimental values are avail- 
able. At high parton momentum transfers, PYTHIA (HO] 
is employed for pQCD scatterings. 

UrQMD differentiates between two regimes for the ex- 
citation and fragmentation of strings. Below a momen- 
tum transfer of Q < 1.5 GeV a maximum of two longitu- 
dinal strings are excited according to the LUND picture, 
at momentum transfers above Q > 1.5 GeV hard interac- 
tions are modelled via PYTHIA. Figure [1] shows the to- 
tal cross-sections of resonant hadronic interactions, string 
excitation and (hard) PYTHIA-scatterings as a function 
of the center-of-mass energy of the collision ^Scoii. The 
contribution of hard scatterings to the total iT~^Tr^ cross- 
section at the highest SPS-energies {^/s ~ 17.3 GeV) is 
about 4 %. Figure[2]shows a comparison between charged 
particle spectra from proton-proton collisions calculated 
in UrQMD with and without the PYTHIA contribution. 
For detailed information on the inclusion of PYTHIA, 
the reader is referred to Section II of [49j . 

In the UrQMD framework, all particle properties 
(mass, width, spectral shape) are taken at their vacuum 
values, the propagation is performed without potentials 
(cascade mode). UrQMD has been used by Dumitru et 
al. to study direct photon emission earlier [46} : a brief 
comparison between their results and the results obtained 
with this approach can be found in Appendix [Cl 
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FIG. 2: (Color Online.) Charged particle spectra from 
proton-proton-collisions at iJiab ~ 158 GeV calculated with 
UrQMD with (solid black line) and without (dashed red line) 
PYTHIA. 



B. Hybrid model 

In the following, we compare results from this micro- 
scopic model to results obtained with a hybrid model de- 
Here, the high-density part of the reaction 
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scription 

is modelled using ideal 3+1-dimensional fluid-dynamics. 
The uncquilibratcd initial state and the low-density final 
state arc described by UrQMD. Thus, those stages are 
mainly governed by string dynamics (initial state) and 
hadronic rescattering (final state). 

To connect the initial transport phase to the fluid- 
dynamic phase, the baryon- number-, energy- and 
momentum-densities are smoothed and put into the hy- 
drodynamic calculation after the incoming nuclei have 
passed through each other. Note that in non-central 
collisions, the spectators are propagated in the cascade. 
Thus, the initial state for the hydrodynamic stage is 
subject to both geometrical and event-by-event fluctu- 
ations. Temperature, chemical potential, pressure and 
other macroscopic quantities are determined from the 
densities by the Equation of State (EoS) used in the cur- 
rent calculation. During this transition, the system is 
forced into an equilibrated state, regardless of the actual 
level of equilibration before the transition. The initial 
temperature profile at z = for a sample Pb + Pr- 
event with & = fm and i?Lab ~ 158 AGeV for two 
different EoS (Bag Model and Hadron Gas) is shown in 
Figure [3l Then the ideal (3-1-1) dimensional hydrody- 
namic equations are solved on a grid using the SHASTA- 
algorithm [36j . 

After the local rest frame energy density has dropped 
below a threshold value of Ccrit = 730 MeV/fm"^ (« 5eo), 
particles are created on an isochronous hyper-surface 
from the densities by means of the Cooper-Frye formula 
and propagation is continued in UrQMD. 

The transition scenario chosen for the present studies 
is always isochronous, i.e. the whole system must meet 
the criterion at the same calculational-frame time before 
the transition is performed. Earlier investigations within 
this hybrid model include an extensive analysis of the ef- 
fect of changing the transition criterion [5l| , strangeness 



-5 - 



-10 



BM-EoS 



1 HG-EoS 




10 12 



FIG. 3: (Color Online). Temperature profiles after switch- 
ing to fluid-dynamic description, calculated with the Bag 
model (left half) and hadron gas Equation of State (right 
half). Along the lines, the temperature is constant, going 
from T = 50 MeV at the outermost line to T = 250 MeV 
at the innermost. Calculations have been done for Pb + Pb- 
coUisions at i?Lab = 158 AGeV with 6 = fm. 



production [5^ |53|, HBT correlations 
mass [5^ and elliptic fiow [56l |. 



transverse 



C. Equations of State 

For the investigations presented here, different Equa- 
tions of State are used for the hydrodynamic phase. The 
base line calculations are done with a hadron gas Equa- 
tion of State (HG-EoS), which includes the same degrees 
of freedom as present in the transport phase. This allows 
to explore the effects due to the change of the dynamic 
description. Secondly, a MIT-Bag model EoS (BM-EoS) 
with a partonic phase and a first order phase transi- 
tion [3^ is employed. The BM-EoS thus allows for in- 
vestigations of photon emission from the QGP. In order 
to obtain meaningful values of temperature and chemical 
potentials from the densities, the BM-EoS is smoothly 
transferred to the HG-EoS just above the transition en- 
ergy density. 



III. PHOTON EMISSION SOURCES 

Photon emission is calculated perturbatively in both 
scenarios, hydrodynamics and transport, because the 
evolution of the underlying event is not altered by the 
emission of photons due to their very small emission 
probability. The channels considered for photon emission 
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may differ between the hybrid approach and the binary 
scattering model. Emission from a Quark-Gluon-Plasma 
can only happen in the hydrodynamic phase, and only if 
the Equation of State used has partonic degrees of free- 
dom. Photons from baryonic interactions are neglected in 
the present calculation. Emission of hard photons from 
early pQCD-scatterings of nuclcons is calculated sepa- 
rately and incoherently added to the simulated spectra. 



A. Photons from microscopic collisions 

In the transport part of the (hybrid) model, each scat- 
tering is examined and the cross-section for photon emis- 
sion is calculated. Here, we employ the well-established 
cross-sections from Kapusta et al. and Xiong et 
al. Kapusta and collaborators based their calcula- 

tions on the photon self-energy derived from a Lagrange 
density involving the pion, p and photon-fields 



1 



1 



(1) 

Here, $ is the pion field, Pfj,^, = df,p^~d^Pfj, and = 
diyAp are the p and photon field-strength tensors 
and Df^ = 9^ — ^eA^ — igpP^ is the covariant derivative. 
The p decay constant gp is calculated from the total width 
FjQj of the p meson: 



dpA 



9l 



487r- 



4to?. 



(2) 



The differential cross-sections used for the present in- 
vestigation [ill, [l^l are given in Appendix |D] 

All scatterings during the transport phase are exam- 
ined in order to obtain direct photon spectra. For every 
scattering that may produce photons (i.e. those that have 
initial states equal to the processes listed in Appendix iPj) . 
the corresponding fraction of a photon. 



f^tot 



(3) 



is produced. Here, crtot is the sum of the total hadronic 
cross-section for a collision with these ingoing particles 
(as provided by UrQMD) and the electromagnetic cross- 
section CTem as calculated by the aforementioned formula. 
In order to obtain the correct angular distribution of the 
produced photons and to enhance statistics, for each scat- 
tering many fractional photons are created that populate 
all kinematically allowed momentum transfers t. In this 
procedure, each photon is given a weight AA^^ according 
to 



AAT* 



dt 



-{s,t)At 



O-tot(s) 



(4) 



and the photons are distributed evenly in the azimuthal 
angle ip. The integral tTem(s) = / '^'^^'^/dtdt is performed 
analytically for each channel. The resulting formulae are 
shown in Appendix [E] 

Since the width of the p-meson is not negligible, its 
mass distribution has to be taken into account. For the 
processes with a p-meson in the initial state, the actual 
mass TUp = ^/pjj^ of the incoming meson is used for 
the calculation of the cross-section. If there is a p-meson 
in the final state, then first the mass of the p is chosen 
randomly according to a Brcit-Wigner distribution with 
mass-dependent width. This mass is then used for all 
further calculations of this process. Figure |4] shows the 
cross-sections of the channels listed above as a function 
of ^/s. 



B. Photons from hydrodynamics 



In the hydrodynamic phase photons arc produced frac- 
tionally from every cell on the hydrodynamic grid whose 
energy density is above a threshold ethr = 10~^^ Eq using 
the parametrizations by Turbide, Rapp and Gale j57l |. 
They use an effective non-linear cr-model Lagrange den- 
sity in which the vector and axial vector fields are im- 
plemented as massive gauge fields of the chiral U{3)l x 
U{i)R symmetry to obtain the rates. For details on this 
ansatz, the reader is referred to the original publication 

As mentioned earlier, the processes calculated by Tur- 
bide et al. differ from those considered by Kapusta et al. 
Only the processes tttt — > jp and np — > 777 are therefore 
common in both models. The rate of Turbide et al. for 
TTp — > 77r directly includes the process with an interme- 
diate ai-meson. 

To simplify the calculations, all photon rates in [s^l 
are parametrized by the general form 



^dR , f B E 



(5) 



where A, C and D are linear functions of some power 
of the temperature T: A{T) = Ax + A2T^\ The pa- 
rameter set can be obtained from 571. In the rates, the 



energy E and temperature T are to be given in units of 
GeV, and the result will have the unit GeV^^ fm^''^. We 
also employ the hadronic form factor introduced in [57| . 

In the Quark-Gluon-Plasma, the rate used is taken 
from Ref. 
result as 



58|. They computed the full leading-order 
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FIG. 4: (Color Online.) Cross-sections for all included channels as a function of ^/s. For visibility, the cross-sections for all 
processes np — > "fn are shown separately. They have been calculated with a p mass mp = 0.769 GeV. The left plot shows the 
cross-sections for tttt — >■ 7p both for fixed p mass {nip = 0.769 GeV, labelled "nip fixed") and for variable p mass (labelled "nip 
Breit-Wigner" ) . 



dR ^ 2 Q:emas ^2 1 A ( ^\ , 1 



In (2a;) + C22(x) + Cb 



(6) 



and give convenient parametrizations for the contribution of 2 <-> 2-, bremsstrahlung- and annihilation-processes (C22, 
Cbrcms and Cann, respectively) 



^brems('^) ^~ C'ann(*^) 



C22{x) = 0.041X-1 - 0.3615 -f 1.01 cxp (-1.35a;) 

lv7 [0.548 In (12.28 + i) 0.133x 



6.27 



(7a) 
(7b) 



In Equations ([S]) and x — ^/t, qi is the charge of 
quark- flavour i, acm and ag = 9^ /in are the electromag- 
netic and QCD coupling constants, respectively. In our 
calculations, we use Nj — 3, and therefore J2i 1i = V^- 
The temperature dependence of as is taken from [59| as 



as(T) 



(33-2iVy)ln(f )' 



(8) 



and the critical temperature at fig = to be Tc 
170 MeV. 



C. Photons from primordial pQCD-scatterings 

At high transverse momenta, a major contribution to 
the photon yield is the emission of photons from hard 
pQCD-scatterings of the partons in the incoming pro- 
tons. In the intermediate and low p_L-regions, the contri- 
bution may be comparable to or smaller than the yield 
from other sources. 

We apply the results extracted by Turbide et al. \E% . 
They first scale the photon spectrum from proton-proton- 
collisions by the number of binary collisions in Pb-|-Pb- 
coUisions, and then add a Gaussian-shaped additional 
fcj^-smearing to the result. The width of the Gaussian 
is obtained by fitting this procedure to the data from 
proton-nucleus collisions. The results shown here are ob- 
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FIG. 5: (Color Online.) Comparisons between the rates from 
[s?! (lines) and box calculations restricted to a tt, p, ai-system 
(points) with UrQMD at T = 150 MeV. 



tained with a (Afc^ ) ^ 0.2 GcV^. 

For comparison, we also show pQCD spectra obtained 
earher by Gale [sO] following Wong et al. [6l|. They 
follow the same procedure as explained above. The au- 
thors of [60} obtain a higher intrinsic transverse parton 
momentum of (Afc^) = 0.9 GcV^, yet lower spectra. 

In order to compare our calculations to experimental 
data in Figure [HI we use the newer calculations by Tur- 
bide et al. 57 1 . 



V. PHOTON SPECTRA 
A. Emission stages 

For the present investigation, the evolution of a heavy- 
ion reaction is divided into three stages: "initial stage" , 
"intermediate stage" and "final stage" . The time spans of 
these stages are defined by the calculations in the hybrid 
model as before {t < 1.4 fm), during (1.4 fm < t < 
13.25 fm for the hadron gas and cascade calculations and 
1.4 fm < t < 30 fm for the Bag model calculations), and 
after the hydrodynamic description (t > 13.25 fm (HG, 
cascade) and t > 30 fm (BM)). I.e., the initial phase 
denotes the initial string dynamics until the transition to 
hydrodynamics. Since all differences between the models 
start at the transition to hydrodynamics, this stage and 
its contribution is always the same. 

Secondly, the "intermediate stage" denotes the phase 
described by hydrodynamics. It starts with the transition 
to hydrodynamics and ends with the transition back to 
the cascade. Within this stage, the degrees of freedom 
may be partonic, hadronic and partonic (in case of a Bag 
model EoS calculation) or hadronic (in all cases). The 
"final stage" starts at the transition from hydrodynamics 
to the cascade. On average this happens at 13.25 fm for 
the hadron gas EoS and 30 fm for the Bag model EoS. For 
pure transport calculations, this phase is also set to start 
at 13.25 fm, so that a comparison between hadron gas 
EoS and cascade calculations is possible. The degrees of 
freedom in the "final stage" are hadrons and resonances. 
The calculations are done with UrQMD in cascade mode, 
and proceed until the last collision has happened. 



IV. RATES FROM TRANSPORT AND 
HYDRODYNAMICS 

Before comparing photon spectra from complex 
nucleus-nucleus collisions between cascade- and hybrid 
model, we check if both approaches give similar results 
for the setup of a fully thermalized box. 

I.e., we perform UrQMD calculations in a box (6^, 
allowing only tt-, p- and ai-mesons to be present and to 
scatter. When the matter in the box has reached thermal 
and chemical equilibrium, the rate of photon emission is 
extracted based on the microscopic scatterings with the 
procedure described in Section fill Al Then, we compare 
the microscopic rates to the hydrodynamic rates from 
Equation ([5]) with the parameters from [53|. Since the 
available rates in the cascade and hydrodynamic modes 
differ, as pointed out above, we restrict the comparison to 
the common rates tttt — >■ 7P and tt/? — > jn. The cascade- 
rates are explicitly summed over all charge combinations. 

Figure [5] shows the comparison of the rates at a tem- 
perature T = 150 MeV. It can be seen that the micro- 
scopically obtained rates agree very well with the ther- 
modynamic rates. 



B. Influence of EoS and dynamics 

The contributions of each of these phases to the fi- 
nal spectra are shown in Figure [6] for the pure cas- 
cade calculation and in Figure [7] for the hybrid calcu- 
lation with the hadron gas EoS. In both figures, the 
black solid lines show the complete direct photon spec- 
tra, the green dash-dotted lines show the initial stage- 
contribution, the red dashed lines show the photon spec- 
trum from the final phase and the blue dotted lines show 
the contribution from the intermediate stage. In the hy- 
brid model, the evolution during the intermediate stage 
is calculated within the hydrodynamic framework. We 
find that the contribution from this intermediate phase 
(1.4 fm < t < 13.25 fm) is similar in both models. How- 
ever, at high transverse momenta > 3 GcV, the cas- 
cade calculation yields more photons from this interme- 
diate time span. This can be related to imperfect ther- 
malization of the system at the transition from the initial 
non-equilibrium state, which is forced to thermalization 
at the transition to the hydrodynamic phase, but pre- 
served when doing cascade-only calculations. However, 
Figure [S] suggests photon emission towards high trans- 
verse momenta from the intermediate stage is in any case 
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FIG. 6: (Color Online.) UrQMD calculation. Contributions 
of the initial {t < 1.4 fm, dash-dotted green line), intermediate 
(1.4 fm < t < 13.25 fm, dotted blue line) and final (13.25 fm < 
t, dashed red line) stages to the spectrum from pure cascade 
calculations (solid black line). 



FIG. 8: (Color Online.) Hybrid model calculation with BM- 
EoS. Contributions to the initial (green dash-dotted line), in- 
termediate (blue dotted line) and final (red dashed line) stages 
to the inclusive spectrum (solid black line). In the intermedi- 
ate stage, the matter is described by hydrodynamics. 
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FIG. 7: (Color Online.) Hybrid model calculation with HG- 
EoS. Contributions of the initial (green dash-dotted line), in- 
termediate (blue dotted line) and final (red dashed line) stages 
to the inclusive spectrum (solid black line). In the intermedi- 
ate stage, the matter is described by hydrodynamics. 



strongly suppressed with respect to photon emission from 
the initial stage. It is therefore justified to neglect non- 
equilibrium effects from the intermediate phase in the hy- 
brid model. The major difference between cascade and 
hybrid model calculations is the magnitude of the con- 
tributions from the final stage (red dashed lines). Here, 
both models describe the system in the same way. The 
contribution from this phase to the hybrid model direct 
photon spectrum is very similar to the contribution from 
the intermediate stage. In the pure cascade calculation, 
on the contrary, the final phase contributes roughly a fac- 
tor of 5 less to the spectrum than the intermediate stage 
does. A possible explanation for this is that the transi- 



tion procedure from hydrodynamic to transport descrip- 
tion enhances the number of meson-meson collisions in 
the late phase. 

Interactions at early times {t < 1.4 fm) are a signifi- 
cant source of high transverse momentum p± > 3 GcV 
photons. Here, the thcrmalizcd system of the later phase 
provides only minor contributions to the inclusive spec- 
trum. 

In the same analysis done for the hybrid model with 
the Bag model EoS (see Fig. [8]), the picture is differ- 
ent. The contribution from the hydrodynamic interme- 
diate stage is strongly enhanced, and the contribution 
of the final stage after the transition from the hydrody- 
namic phase is, in turn, reduced with respect to the cases 
presented before. The total photon spectrum is com- 
pletely dominated by emission from the hydrodynamic 
phase at low and intermediate transverse photon mo- 
menta p± < 3 GeV. This coincides with the observation 
that the length of the hydrodynamic phase in this model 
(sa 29 fm) is much longer than in the calculations with 
hadron gas EoS (w 12 fm), due to the large latent heat 
in the setup with the Bag model EoS. At high transverse 
photon momenta, the initial stage non-equilibrium cas- 
cade phase dominates the spectrum, as in the other cases 
discussed above. 



C. Channel decomposition 

The contributions of the different channels to the 
hadronic calculations, both in the pure cascade mode 
and the hybrid approach with a hadron gas EoS, is 
shown in Figures [5] and 1101 The dominant contributions 
in both models are very similar; at low transverse mo- 
menta p± < 0.5 GeV, the process with two initial pi- 
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FIG. 9: (Color Online.) UrQMD calculation. Contributions 
of the different channels. 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 

P± [GeV] 



FIG. 10: (Color Online.) Hybrid model calculation with 
hadron gas EoS. Contributions of the different channels. 



ons TTTT 7p is dominant, while in the broad range of 
0.5 GcV < < 3.5 GeV, irp — >■ jir is the major source 
of photons. Processes with an eta-meson (tttt — > 777 and 
TTT] — > 77r) become important at high transverse mo- 
menta, where they contribute in similar magnitude as 
the two aforementioned channels. Processes with kaons 
contribute less to the spectrum than the corresponding 
non-strange channels, and the process tttt — > 77 is the 
least significant in all calculations, as expected. 

In the photon spectrum extracted from the cascade 
calculation, one can observe a flattening of the spectrum 
at high transverse momenta p± ^ 3 GeV. At this point, 
TTTT — >■ 7/9 and 77-processes start to provide significant con- 
tributions to the photon spectrum and lead to a flatter 
slope already at pj_ w 2 GeV. In Section IVI[ we come 
back to this slope change and show it to be consistent 
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FIG. 11: (Color Online.) Hybrid model calculation with Bag 
model EoS. Contributions of the different channels in hybrid 
calculation with Bag model EoS. 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 

P± [GeV] 



FIG. 12: (Color Online.) Comparison of direct photon spectra 
from the hybrid model with hadron gas EoS and transport 
model, using only common channels. 

with the average emission times of the photons at these 
transverse momenta. 

The calculations with the Bag model EoS yield a differ- 
ent picture (see Figure [TT|) . Here, the dominant contribu- 
tion comes from the Quark Gluon Plasma, whose emis- 
sion magnitude is about two times higher than the com- 
bined contribution from all hadronic processes. Again, 
initial stage (pre-equilibrium) processes are dominant at 
high transverse momenta p_L > 3 GeV. Apart from that, 
we observe a similar distribution among the reaction 
channels as we did above (sec discussion of Figure [T0|) : 
the process irp — > 777 is dominant in the intermediate pj^- 
region, tttt — jp dominates the hadronic contribution at 
low p±, the early scatterings in those channels and the 
channels involving Ty-mesons contribute to the high p±- 
region in approximately equal amounts. 

Figure [T^] shows a comparison of hybrid model calcu- 
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FIG. 13: (Color Online.) Comparison of the direct photon 
spectra from all variations of the model to the experimental 
data by the WA98-collaboration @|. Calculations without 
intermediate hydrodynamic stage (pure cascade mode) are 
shown as red crosses, hybrid calculations with hadron gas 
EoS as red solid line and Bag model calculations are depicted 
by the dark-green dotted line. The data contain no photons 
from initial proton-proton collisions. 



lations with hadron gas EoS and transport-calculations, 
using only common channels applied to a dynamic system 
(see Figure [S] for a static comparison). One observes that 
the direct photon spectrum is not sensitive to the change 
in the underlying dynamics (e.g. finite viscosities vs. ideal 
fluid) and indicates that even the two-body collision dy- 
namic in UrQMD drives the system into equilibrium in 
the TT — p — fli channel. 



D. Comparison to data 

A comprehensive comparison between the models pre- 
sented above and experimental data is shown in Fig- 
ure [131 Here, we show the sum over all channels and all 
stages in the cases of pure UrQMD-calculations and hy- 
brid model with hadron gas and Bag model EoS. The Bag 
model Equation of State yields the highest photon spec- 
tra and provides a reasonable description of the data from 
the WA98-collaboration Q. The purely hadronic cal- 
culations, with and without intermediate hydrodynamic 
calculation, give smaller yields than the scenario with a 
QGP. 

After adding the pQCD-spectra as extracted by [57| . 
we obtain Figure [TH Due to the rather large pQCD- 
contribution the difference between the spectra with 
varying EoS is reduced, and all calculations agree very 
well with the data. 
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FIG. 14: (Color Online.) Comparison of the direct photon 
spectra from the model, with added pQCD-photons from [stI ]. 
to the experimental data by the WA98-collaboration . Cal- 
culations without intermediate hydrodynamic stage (pure cas- 
cade mode) are shown as red crosses, hybrid calculations with 
hadron gas EoS as red solid line and Bag model calculations 
are depicted by the dark-green dotted line. 
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FIG. 15: (Color Onhne.) UrQMD calculation. The px- 
spectra for all photons (solid black line), photons from 
boosted collisions (dotted blue line) and collisions with high 
center-of-mass energy (dashed red line). The grey-shaded 
area indicates the range of the px-region used for the solid 
black lines in Figures [TBI and [T71 



VI. SENSITIVITY TO DIFFERENT EMISSION 
TIMES AND PROCESSES 

To investigate the major sources for photons to the 
transverse momentum spectrum and to explore the sen- 
sitivity to the different reaction stages, we investigate the 
origin of the change of slope of the spectra at high trans- 
verse momenta sa 3 GeV in Figures [T21 [H] and [T7] 
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FIG. 16: (Color Online.) UrQMD calculation. The num- 
ber of photons as a function of the transverse center-of-mass 
momentum of the elementary collision g'^""'. We show all 
photons (blue dotted line), photons that have high transverse 
momentum (black solid line) and photons from collisions with 
high center-of-mass energy (red dashed line). The light-blue- 
shaded area indicates the range of the -region used for the 
dotted blue lines in Figures [T51 and [T71 
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FIG. 17: (Color Online.) UrQMD calculation. The num- 
ber of photons as a function of the center-of-mass energy of 
the elementary collision. We show all photons (red dashed 
line), photons that have high transverse momentum (black 
solid line) and photons from collisions of with high transverse 
center-of-mass momentum. The light-red-shaded area shows 
the range of the -^s-region used for the red dashed lines in 
Figures [15] and 1161 



within UrQMD. Two processes may contribute to high 
p^-photons: (I) collisions with a large ^/s in elementary 
reactions (i.e. early stage collisions) and (II) collisions of 
particles with large transverse flow q±^ but rather small 
•\/s (i.e., late stage collisions). To disentangle these ef- 
fects, we determine the contributions of scatterings with 
high center-of-mass energy ^/s > 4 GeV and high center- 
of-mass transverse momentum gj. > 3 GeV. Figure 1151 
shows the transverse momentum spectrum of photons 
split up into the two contributions. Nearly all photons at 
high transverse momenta (grey-shaded area) come from 
collisions with high center-of-mass energies, whereas the 
contribution of high center-of-mass transverse momenta 
only shows a trivial structure at p± sa gthroshoid _ 3 Q^y 

One cross-check for Figure [TS] is shown in Figure [T51 
Here, we observe that the collision spectrum is exponen- 
tial. Only at high transverse boosts, a deviation from 
an exponential spectrum can be seen. This indicates 
that most photons with high transverse momentum come 
from unboostcd collisions with q± < 1 GeV. The center- 
of-mass energy and transverse momentum of a collision 
show no correlation. 

The latter can be seen in Figure [T71 where the photon 
production rate is shown as a function of the center-of- 
mass energy of the individual collisions. The figure con- 
firms the notion that most photons with high transverse 
momenta come from collisions with high center-of-mass 
energies. It is worthwhile to observe, however, that start- 
ing at ^Scoii. = 7 GeV, each elementary collision pro- 
duces essentially only photons with transverse momenta 
p_L > 3 GeV. 
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FIG. 18: (Color Online.) UrQMD calculation. Average pho- 
ton emission times as a function of the transverse momentum 
for all photons, np — >■ yn-, tttt — >■ 7p- and other processes 
(black solid, red dashed, blue dotted and green dash-dotted 
lines, respectively). The shaded areas correspond to the p±- 
regions used for the curves shown in Figure [T^ 



The distribution of center-of-mass energies also shows 
that the vast majority of collisions happen around the p- 
and a 1 -pole masses. 

Figure [m shows the average emission times of photons 
as a function of the transverse momentum for the var- 
ious channels. It is noticeable that over a very broad 
momentum range, 0.3 GeV < p_l < 2.1 GeV, the av- 
erage emission time stays at a constant level of about 
(iomission) ~ 8 fm. This coincides with the region where 
the process irp — > dominates. Only at high transverse 
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FIG. 19: (Color Online.) UrQMD calculation. Emission time 
distribution of photons from yrp-scatterings for the different 
photon transverse momentum regions indicated in Figure [181 
The vertical lines indicate the average emission time in the 
corresponding px-bin. 



momenta, the early times dominate. This is consistent 
with the findings in Section |Vl that the spectrum clearly 
shows two different temperatures, one in the region be- 
low p±^ = 2.5 GeV and a different temperature above 
> 3 GeV. It also explains why the spectral contribu- 
tions from TTTT — > 7/9 and 77-processes show a much flatter 
slope already &t p± k, 2 GeV. At late times, when the 
average center-of-mass energy has decreased, photons are 
predominantly produced at low transverse momenta. 

In Figurc[T9l we show the emission time distribution of 
photons in various p^-bins for the -np — >■ 77r-processes. It 
is interesting to see that at all transverse momenta, there 
is an initial flash of photons emitted at very early times 
t w 1 fm. In the low pj^-bin 1 GeV < p± < 1.5 GeV, 
a very strong contribution from the bulk emission in the 
hot and dense stage between t = 4 fm and t = 12 im 
raises the average emission time. In the intermediate 
p_L-region 2 GeV < p± < 2.5 GeV, the bulk contribution 
is greatly reduced and shines less bright than the initial 
flash. In the highest pj^-region 3 GcV < p± < 3.5 GeV, 
the late bulk contribution is small and the initial stage 
dominates. However, one should note that due to the 
long lifetime of the intermediate stage, the average emis- 
sion times are shifted to higher values. 



VII. SUMMARY 

In this work, we have studied direct photon emission 
from hadronic and partonic sources within three dif- 
ferent dynamical models. In Section IIIIl we presented 
our model for photon emission from microscopic colli- 
sions f Subsection nil A} and from thermal rates (Subsec- 



tion |TTTB|. Then, we introduced the cross-sections and 
thermal rates used for the present calculations. 

In Section llVi we showed that the emission rates from 
a thermalized microscopic cascade calculation agree very 
well with the thermal rates used in the hydrodynamic 
part of the present model. We discussed the contribu- 
tions from the different stages before, during and after 
the high-density part of the evolution to the direct pho- 
ton spectra for cascade calculations as well as hadron 
gas and Bag model hybrid calculations. It was found 
that the relative contributions of photons in the hybrid 
calculation with hadron gas EoS and cascade simulations 
are similar. In contrast, in the Bag model calculations 
the intermediate high-density hydrodynamic phase takes 
a much longer time and contributes substantially more 
to the photon spectra. 

Investigations that differentiate between the different 
channels showed that the process irp — > jn is the domi- 
nant hadronic source for photon production at intermedi- 
ate photon transverse momenta 0.5 GeV < p± < 3 GeV, 
while TTTT — > 7p is dominant at low photon transverse 
momenta p± < 0.5 GeV, only. 

When comparing the different variations of the model 
(see Fig. [T31), we found that both the hybrid model and 
the cascade model can explain the spectra measured 
by the WA98-collaboration if pQCD-photons are in- 
cluded. By comparison of Figures [T2l and [T3l we can also 
conclude that the photon yields from the hybrid-model 
calculation with hadron gas Equation of State and pure 
transport calculation are equal within uncertainties, if 
the same sets of channels is used in the calculations. 

We also found that photons at high transverse mo- 
menta p± > 3 GeV show a significantly flatter slope 
(and therefore higher effective temperature) than pho- 
tons with lower p±. This effect was attributed to higher 
center of mass-energies that produce these photons, and 
we find that it is not significantly infiuenced by elemen- 
tary collisions that have a high center of mass-boost in 
transverse direction. The analysis of average photon 
emission times showed that photons at pj_ > 3 GeV are 
emitted significantly earlier than at lower transverse mo- 
menta. 

We also discussed different model assumptions (see Ap- 
pendix [X] and |B]), namely photon emission from collid- 
ing string ends and compared those to pQCD-spcctra. 
The result indicated that direct photon emission from 
string ends is restricted to the early phase of the collision 
where the medium is hot. Neglecting those collisions low- 
ers the effective temperature at high transverse momenta 
p± > 2.5 GeV. In any case, string end contributions can 
be neglected in comparison to the pQCD contribution. 

We found that the inclusion of a finite width of the 
produced p-mesons is relevant only at low transverse mo- 
menta. I.e., only in the lowest p_L-bin, the effect of assign- 
ing the pole mass or a mass chosen randomly according 
to the spectral function is important. 
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VIII. OUTLOOK 

The good agreement between the calculations pre- 
sented here and the experimental data shows that the 
cascade+hydrodynamic hybrid model provides an excel- 
lent tool to explore the properties of QCD-matter at en- 
ergies where the (onset of) dcconfincment is expected. 
We plan to extend the investigation to a more realistic 
equation of state with chiral restoration, critical endpoint 
and rapid cross-over [6^. Also, the decays of short-lived 
mesons and baryons which cannot be subtracted by the 
experiment may play a major role in enhancing the direct 
photon spectra and will therefore be investigated. 

Furthermore, direct photon emission will be investi- 
gated for the systems studied or planned to be studied 
at RHIC (PHENIX), LHC (ALICE), FAIR (CBM) and 
SIS-100 (HADES). 

More differential observables like direct photon elliptic 
flow can also be explored in the present model. 
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FIG. 20: (Color Online.) UrQMD calculation. Total yields 
of photons with (dashed line) and without (solid line) the 
contributions from colliding string ends. For comparison, the 
spectrum obtained from PYTHIA is shown as well. 
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nential and docs not exhibit the flattening at high trans- 
verse photon momenta. The inclusion of (di-)quark scat- 
terings, however, leads to a strong increase of the photon 
yield at high p±. The contribution of pQCD-photons to 
the inclusive spectrum starts to be significant already at 
relatively low transverse photon momenta p± fa 1 GeV, 
although the magnitude of the contribution differs be- 
tween the different parametrizations in 57 1 and [60j . 



Appendix A: Discussion: String ends 

In UrQMD, the leading particles from a string have a 
reduced cross-section during their formation time. For all 
other purposes, they are treated as hadrons. (Hadronic) 
scatterings of the string ends happen typically at high 
center-of-mass energies ^/s. Collisions from string ends 
can only produce photons if the collision of fully formed 
hadrons of the same type would produce photons. Thus, 
their contribution is treated as an addition to the prompt 
contribution from primordial nucleus-nucleus interac- 
tions. Photons from those collisions contribute signifi- 
cantly to the spectra at high transverse momenta. Since 
these particles arc not fully formed hadrons, but effec- 
tively represent quarks or di-quarks, a hadronic treat- 
ment of those processes is questionable. 

The effects of including the photons from colliding 
string ends, i.e. interaction of leading (di-)quarks, in the 
calculation can be seen in Figure 1201 The spectrum ob- 
tained by neglecting the collision of string ends is expo- 



Appendix B: Discussion: p-meson width 

Earlier, we discussed the handling of the p-meson's 
finite width. Figure [21] shows the effects of follow- 
ing the calculation outlined there. In both channels, 
tt'^tt^ — > 7p° and tt'^tt^ — >■ 7/9^, the yield is about 10 % 
higher for ps produced at their pole mass, and only at 
very low momenta this excess becomes as large as 40 %. 

This behaviour can be explained by kinematic argu- 
ments: The by far highest scattering rate in n + n- 
coUisions is at y/s « m^. Here, the photon cross-section 
with fixed pole mass is much higher than the extended 
calculation with variable p mass. At all other center-of- 
mass energies, the extended model gives a higher cross- 
section, but these comparatively rare processes provide 
only a minor contribution to the spectrum. The processes 
at low will contribute primarily to the low-p_L-region, 
because the production of the p-meson consumes most of 
the available energy. Therefore, the enhancement in the 
model is most pronounced at low p±. 
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FIG. 21: (Color Online.) UrQMD calculation. Photon spec- 
tra from collisions that produce a p meson in the final state 
for the production of p mesons at its pole-mass (solid and 
dashed lines) and for the production of p mesons according 
to a Breit-Wigner-mass distribution (dotted and dash-dotted 
lines) . 
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FIG. 22: (Color Online.) Comparison of the current model 
(red dashed line) to calculations from Dumitru et al. ([4^. 
black solid line) . For further comparision, we apply our model 
to an earlier UrQMD version (blue dotted line). 



Appendix C: Comparison to older works with 
UrQMD 



As mentioned in Section|TTl UrQMD has been used pre- 
viously in order to obtain direct photon spectra . The 
authors calculated transverse momentum spectra for cen- 
tral (6 = fm) Pb-fPb-colhsions at £:Lab = 160 AGeV. 
Limiting ourselves to the same conditions, wc can com- 
pare our work to that of Dumitru et. al. In Figure !^ we 
compare results from the current UrQMD- version 2.3 to 
those obtained with the earlier versions 1.3 (using our 
own photon-analysis) and 1.0 (taken from [4^). The 
older UrQMD-versions yield significantly higher photon 
spectra at intermediate p_L . 



We can identify two changes in UrQMD that lead to 
this behaviour. Between versions 1.0 (used by Dumitru 
et al.) and 1.3, the angular distributions of various pro- 
cesses have been altered. Since the collisions of pions 
with a high difference in rapidity provide significant con- 
tributions to the spectra from [4y| , the improved angular 
distributions reduce the photon production cross-section. 
Furthermore, in versions prior to 2.3, the number of pi- 
ons has been unphysically high. The correction of this 
leads to fewer collisions involving pions and hence to a 
further reduction of the spectra. 



Appendix D: Differential cross-sections 



For completeness, we list the differential cross-sections 
for the processes used in [ll|, [l3] ■ 
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In these equations, t = {p-^ — p^)^ is always the 
momentum transfer from the pion to the photon [64 1 
(unlike the convention used in [Tl| ) and Pc.m. = 
(2^5)" ~ 2s(777^ + 777^) + (777^ - 777^)2 is the three- 
momentum of the incoming particles in the center-of- 
mass frame (7711 and 7772 being their masses). The value 
of A is, consistent with [ll|, A = g^ppgp/4:TTjp = 4.7, and 
the pion electromagnetic form factor is 



In their 1992 paper, Xiong et al. (l2[ calculate the 
cross-section for the formation of an intermediate oi- 
meson during Trp-scattering, averaged over all possible 
charge combinations: 
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This channel is not included in Kapusta et al. Xiong et the pion-, photon-, p- and ai-fields 
al. obtain this from a Lagrangc-density involving only 
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C Gpa^{gp^{p^ -pp) - p-KpPp^,)p"^ + Gp—a^{gf,^{prr ■ Pj) ~ p^pPj^)A''<^ , (D2) 
where Gp = 14.8 GeV^^. The partial widths of the ai, Ta-^^Trp and Ta-i^^^T^, are estimated to be 
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Appendix E: Integrated Cross-sections 

The integration over the cross-sections hsted in ap- 
pendix |D] yield the fohowing resuhs: 
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, (Eld) 



ml - -sy^itl tl) ^-{tl ti) ml {mt^ + s{ml - m^)) 



{2ml 

nl [{t- + u^){s - ml) + {2ml ~ s)^] h 

(i- + - ml)Io + hi 



(Ele) 

m) 



S [At + (i_ + 7i_ - 7772)/i + {t_ + 7i_ - ?772)2/o - 7772r2/o] 



27ra2 

SPcm. 

Ap,^- {np 



At + 2777, 



fll 



1 

77r) 



27772 



In- 



4- I 

t+ - ml 



m-iAt 



(t- -ml){t+ - ml) 



+ {t± ^ m=f) 



(Elg) 
{^Ih) 
(Eli) 



In those equations, t± = m\ — 2iu{Ei ± pi) are the correspondingly u± = rnl + m^ + m^ — s — t±, and At is a 
minimal and maximal allowed momentum transfers, and 
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shorthand for i_ — t+. The indices 1, 2 and 4 denote the 
particles in the order given on the left hand side of the 
equations (3 is the photon). Ei, pi and u ^ denote 
momentum and energy of the respective particles in the 



center-of-mass frame. In the cross-section tt^t] — ^tt^ 
(Eq. ( |Elg[ )), the following notations have been used for 
simplicity: 



at an 



atan 



and 



h = In 



p p 



(E2a) 
(E2b) 
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